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ROYAL  AIRCRAFT  ESTABLISHMENT ,  FARNBOROUGH 

A  Study  of  Azimuth  Guidance  and  Control 
for  the  Vernier  Stage  of  Blue  Streak 


by 

R.  H.  Merson,  B.A. 
and 

J.  M.  Howl,  B.A. 


SUMMARY 


A  control  system  in  azimuth  based  on  a  small  swivelling  jet  is  proposed 
for  the  vernier  stage  of  the  medium  range  ballistic  missile  Blue  Streak. 

This  control  system  is  applicable  to  both  inertia  and  radar  guidance  systems, 
the  measured  error  quantity  being  essentially  the  azimuth  component  of  the 
deviation  of  the  missile  flight  path  from  the  line  joining  missile  and  target 

The  jet  has  to  be  sufficiently  small  to  avoid  a  missile  speed  error  at 
cut-off.  This  leads  to  a  relatively  small  available  deflecting  foroe  in 
azimuth  and,  to  make  full  use  of  this,  a  system  with  a  high  loop  gain  and 
limited  error  input  is  proposed,  stabilization  being  obtained  with  heading 
and  heading  rate  feedback.  The  input  limits  and  gain3  are  chosen  so  that, 
in  effect,  the  missile  heading  is  limited  within  prescribed  bounds. 

The  performance  of  the  system  has  been  studied  on  an  analogue  simulator 
under  assumptions  of  constant  missile  parameters  and  the  visual  small  angle 
approximations.  It  i3  shown  that,  with  proper  choice  of  oontrol  coefficients 
the  system  is  very  little  worse  than  an  ideal,  infinitely  stiff,  system  with 
the  same  heading  limitation. 

The  steady  state  error  due  to  heading  gyro  wander  is  considered,  and  for 
the  missile  parameters  used  is  thought  to  be  tolerable. 

The  effects  of  the  data  sampling  and  computer  time  lag  relevant  to  the 
radar  guidance  system  have  been  studied  and  it  is  shown  that  sampling  periods 
and  time  lags  up  to  about  1  seoord  have  no  material  effeot  on  the  oontrol 
system. 
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1  Introduction 

The  controlled  flight  of  the  medium  range  ballistic  missile,  Blue 
Streak,  can  be  divided  into  two  propulsion  stages.  During  the  first  or 
main  propulsion  stage,  the  missile  is  boosted  to  a  velooity  of  about 
16,000  ft/seo  at  a  height  of  some  100  h.  miles.  The  overall  noise  in  the 
control  system  during  this  stage  may  be  too  large  for  (a)  sufficiently 
accurate  azimuth  control,  and  (b)  sufficiently  accurate  cut-off  timing. 
Provision  is  therefore  being  made  for  a  seoond  propulsion  stage,  the  vernier 
stage,  using  small  auxiliary  motors.  During  this  stage  the  noise  level  will 
be  lower  and  the  problems  of  removing  the  residual  azimuth  errors  and 
effecting  the  correct  "cut-off"  are  simplified. 

Control  of  the  missile  throughout  the  propulsion  stages  is  obtained 
solely  by  swivelling  the  jets.  During  the  main  stage  the  control  problem 
is  complicated  by  the  rapidly  varying  conditions  of  speed,  mass,  and  aero¬ 
dynamic  forces  on  the  missile.  During  the  vernier  stage,  however,  which  is 
virtually  in  vaouo,  the  aerodynamic  foroes  are  negligible  and  the  variations 
in  speed  and  mass  of  the  missile  are  very  small.  Par  this  reason,  the 
problem  of  azimuth  control  during  the  vernier  stage  has  been  tackled  first, 
and  this  note  deals  with  the  design  aid  simulator  study  of  a  suitable  con¬ 
trol  system. 

As  will  be  shown,  this  study  is  relevant  to  both  inertia  and  radar 
guidance  systems.  In  the  case  of  the  radar  system,  further  work  is  inoluded 
on  the  effect  of  computer  lags  and  the  intermittent  nature  of  the  error 
signal. 

2  Equations  of  motion 

In  Pig.1  Ox  and  Oy  are  taken  as  reference  space  axes  in  the  azimuth 
plane.  The  axis  of  the  missile  makes  an  angle  ijr  (heading  angle)  with  Ox, 
and  the  thrust,  P,  acting  at  a  distance  l  to  the  rear  of  the  centre  of 
gravity,  makes  an  angle  £  with  the  missile  axis.  The  tangent  to  the 
flight  path  makes  an  angle  ijfo  with  Ox  and  the  target  direction  makes 
an  angle  ijrp  with  Ox.  Let  m  be  the  mass  of  the  missile,  B  its  lateral 
moment  of  inertia  and  v  the  magnitude  of  its  velocity  component  in  the 
azimuth  plane. 


The  equations  of  motion  then  are:- 

mX  =  P  oos  (\[h-£)  (1 ) 

my  =  P  sin  (2) 

B$  =  -  P4  sin  £  (3) 

*  =  V  cos  Ijrm  (4) 

y  a  -  v  sin  Vm  .  (5) 


These  can  be  reduced  by  eliminating  x  aid  y.  Differentiating 
equations  (4)  and  (5)  and  substituting  for  5,  $,  in  equations  (l)  aid  ( 2) :  — 


m  (*  cos  sin  i|rj 

=  P  COS 

(6) 

m  (v  sin  \|rm+  v^  cos  hj 

a  P  Sin  (*£) 

(7) 

4  ** 
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whence 

V  =  f;  COS  (**-&-ijrm) 

(8) 

and 

(9) 

and  from  equation  (3), 

i  *  -f  sln^.  ‘  (10) 

Equations  (8)  to  (10)  describe  the  motion  of  the  missile  in  response  to  a 
given  Jet  defleotion  £.  However,  they  are  unnecessarily  complicated  for  the 
order  of  investigation  initially  required.  The  following  assumptions  are 
therefore  made:- 

(i)  £  and.  are  sufficiently  small  for  siii  £,  sin  (^n-£~^m)  to 

be  replaced  by  £, 

P  p/, 

(ii)  —  and  -=r  are  constant.  . 

'  '  mv  B 

P  0 

(From  equation  (8)  and  the  small  angle  approximation,  —  £  v,  so  that  the 
former  assumption  is  equivalent  to  ^  =  constant.) 

With  these  assumptions,  the  equations  of  motion  reduce  to:« 

(11) 
(12) 


For  the  present  study  the  values  T^  a  2000  secs,  wQ  =  1  rad/sec,  were  taken. 

These  could  represent,  for  example,  a  missile  with  a  thrust-to-mass  ratio  of 
:£g,  a  speed  of  16,000  ft/sec,  and  a  moment  of  inertia  (taking  B  *  J  mi2) 
corresponding  to  a  uniform  thin  rod  of  half-length  i  a  24  ft. 


Guidance  and  oontrol  equations 


In  the  case  of  the  inertia  guidance  system  and  one  possible  method  of 
radar  guidance  computation,  the  measured  error  quantity  is  the  line  error,  AL, 
at  the  target,  assuming  the  missile  were  to  continue  along  its  present  flight 
path  (Pig, 2). 
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Thus, 

AL  *  R 

•where  R  is  the  range  to  go. 

Now  R  changes  by 'perhaps  100  miles  during  the  vernier  stage,  that  is 
to  say  by  perhaps  %  of  its  value  at  cut-off.  For  the  present  purpose  it  can 
therefore  be  considered  constant,  and  .  is  proportional  to  the  line 

error  AL. 

In  the  case  of  another  proposed  method  of  radar  guidance  computation, 
the  azimuth  error  is  measured  as  the  lateral  velocity  component  perpendicular 
to  the  desired  velocity  direction,  that  is  to  say,  v  sin  ('Jfm~'|»p5»  But, 
since  v  may  be  considered  constant  and  sin  (tm"*'l,T)  ^  it  follows 

that  is  again  an  adequate  representation  of  the  error  measurement. 

Thus  the  two  methods  of  guidance  in  azimuth  are  identical;  in  the 
inertia  case,  however,  the  estimation  of  the  line  error  AL  is  possible 
only  if  the  missile  is  close  to  a  predetermined  trajectory,  while  in  the 
radar  system  there  is  no  such  restriction. 

In  practice,  the  measured  value  of  (t-^iji)  will  lag  the  true  value, 
either  accidentally,  because  of  instrument  lags,  or  purposefully,  to  intro¬ 
duce  smoothing. 

Let  ijfg  be  the  output  of  the  error  measuring  system,  i.e.  the  guidance 
system,  and  assume  that  this  system  has  a  linear  lag  with  time  constant  T0. 

The  guidance  equation  then  takes  the  form:- 


(13) 


In  the  case  of  a  ground  computer  for  measuring  this  is  not  a 

realistic  approach,  but  as  already  mentioned,  the  effect  of  finite  time  delays 
is  considered  in  a  later  section. 

The  proposed  control  equation,  in  general  terms,  is:- 


%  =  "  KiY1  (D)  .  +  Kgt  + 


(14) 


in  the  usual  notation. 

Here  \jr  and  its  derivative,  are  assumed  to  be  measured  by 
heading  and/or  rate  gyros  in  the  missile.  The  bandwidth  of  the 
motor-swivelling  servos  is  assumed  to  extend  beyond  the  highest  frequencies 
likely  to  be  encountered  in  the  missile  behaviour.  This  is  reasonable  since 
the  system  frequencies  are  relatively  low. 

A  simple  block  diagram  of  the  complete  system  is  shown  in  Pig.  3,  and  a 
representation  of  the  linearized  equations  in  Pig.4« 
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4  Limit  on  missile  heading  angle 

In  this  note  the  term  "heading  angle"  is  occasionally  used  in  contexts 
where  it  might  be  thought  preferable  to  use  the  term  directional  error,  the 
latter  term  being  defined,  in  this  oontext,  as  (^-ta).  However,  if  it  is 
assumed  that  the  heading,  gyro  datum  is  parallel  to  the  target  direction 
(within  a  few  degrees)  then  the  heading  angle,  f,  as  measured  by  the  gyro,  is  a 
sufficiently  good  estimate  of  the  directional  error,  In  a  radar 

guided  system  it  is  desirable  to  limit  the  directional  error  so  that  its  magni¬ 
tude  does  not  exceed  say  30°.  However,  in  the  simulation  of  suoh  a  system, 
with  typ  a  0,  it  is  the  heading  angle  j  ljr,  which  is  limited,  so  that  this  is 
the • situation  treated  in  the  analysis. 

5  Preliminary  analysis  -  the  ideal  system 

.The  main  object  of  the  vernier  stage  (apart  from  cut-off  adjustment)  is  to 
eliminate  azimuthal  errors  as  quickly  as  possible,  at  the  same  time  avoiding 
large  heading  and  jet  angles.  This  means  that  the  oontrol  system  will  be 
highly  non-linear,  and  a  linear  analysis  is  not  sufficient. 

Some  insight  into  the  problem  may  be  obtained  by  considering  a  system 
which  we  may  call  the  ideal  one.  Let  us  suppose  that  ljr  is  limited  in  magni¬ 
tude  to  the  value  '^nax.  Then  in  order  to  obtain  the  best  response  to  a  step 
in  ij»p  (which  is  equivalent  to  having  an  initial  lateral  velooity  error),  the 
best  that  oan  be  done  is  to  put  on  the  maximum  heading  angle  as  rapidly  as 
possible  in  order  to  obtain  the  maximum  lateral  acceleration,  hold  this  heading 
until  the  lateral  velocity  is  correct  (i.e,  ijre  a  0),  and  then  reduce  the  head¬ 
ing.  rapidly  to  a  small,  value.  This  is  illustrated  in  Pig. 5* 

Since  ijr  starts  ideally  with  a  step,  £  (proportional  to  ijr)  is  a 
second-order  6-function.  The  angular  error,  ijre,  ideally  decreases  linearly 
to  zero,  for  when  i(r  a  £  a  0,  ^  «  ljr  and  therefore  $  VT1  •  W 

from  equation  ('ll);  hence,  from  equation  (13),  since  is  oonstant, 

T  **  T,|  ^max  “  **  mv  ^max  *  0.5) 

The  rate  of  reduction  of  error  is  thus  approximately  proportional  to  the 
maximum  permissible  heading— to-f light  path  angle,  and  to  the  thrust- to-mass 
ratio. 

It  will  be  seen  later  that  equation  (15)  is  very  nearly  realised  in 
practice. 

6  The  control  equation 

In  order  to  achieve  a  system  which  approaches  the  ideal  one,  we  are  bed 
to  a  oontrol  equation  involving  heading  and  rate  of  ohange  of  heading, 
together  with  a  very  high  loop  gain  in  the  guidance  loop  aid.  a  limit  on  the 
guidance  input. 


Thus:— 

■  -^(d)  (i&o 

for  |T1(D)l£|<A 

%  »  +K1A  +  K2t+K3t  (l6b) 

for  ^(D)  >  A  . 
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<33  (the  demanded  jet  angle)  is  the  input  to  the  jet  servos;  if  it  is 
assumed  that  the  achieved  jet  angle  is  limited  by  mechanical  stops  to  +£t,, 
then 


«  =  *D  for  M  <^L 

S  »  for  |^|  >  ^ 


neglecting  the  servo  lags,  as  stated  previously. 

7  Preliminary  choice  of  parameters 

The  stability  of  the  system  must  be  considered  in  the  two  separate  cases 
of  settling  to  the  maximum  heading  angle  (the  'non-linear1  settling),  and  the 
final  linear  settling. 

Considering  the  settling  to  the  maximum  heading  angle,  when  the  jet 
deflection  is  off  the  limits  but  (D)  >  A,  we  have 

“  —^2  =  +K,A  +  K^+  Kjf  (18) 

wo 

which  is  a  linear  differential  equation  in  ijr,  with  the  steady  state  heading 
K.A  j. 

angle,  ^  =  +  -g —  ,  undamped  angular  frequency  w  a  wq  K^2,  and  damping 

.  .  2*2*0 
ratio,  u  => - X"  • 

Vs 


Now  we  require  a  steady  state  I*  I  a  *  and  a  damping  ratio  of  0.8, 

4  s 1  nicLX 

say,  which  gives  the  relations 


K.A 

Y~~  =»  ijr 

Kj  Tm 


1.6  K, 


Consider  now  the  final  stage,  taking  (D)  s'  1, 


The  equations  are 


i  =  x 
1 


♦  = 


«  =  (*r*m)  +  \\  + 


Eliminating  £  and  ijrm  one  obtains 
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(D  *  *  »oV +  »oS>  *  ^  -  d2>]  *  -  v.2  (D  *  i-)  % 


where  D  2 


at  * 


The  stability  polynomial  is  thus  a  cubic:- 

2 


■Kjw 


T. 


1 


To 


\  /K,w  2  K.w 


,  2  K2w  2  K,Wq2 

Since  »  1,  T^  »  1 ,  the  terms  -j-  p  ,  — -  p,  ana  -  , 

be  neglected  for  the  present  purpose,  leaving  the  sampler  form:-* 

K.w  2 
1  o 


3/2  K1\  2  2 

^foV^+K3V  +  i; 


can 


(22) 


1.6w 


Finally,  making  the  substitutions  K,  =  — g  ,  Kg  =  -  *-g-  (i.e.  taking  the 

w  w 

o  o 

damping  ratio  of  the  non-linear  settling  to  be  0.8,  as  above),  the  stability 
cubic  becomes 


3  (.  ,  Kl\  2  ■  2_  K1Wo 

p  +  ( 1 . 6w  ■  -  ■ 


-*) 


p  +  w  p  + 


T. 


1 


(23) 


This  cubic  can  be  written  in  the  form:- 

(p2  +  2<xp  +  P)  (p  +  r) 

and  the  roots  were  first  investigated  for  the  oase  wQ 
sec,  T^  b  2000  secs  and  various  values  of  ;  the 
Figs.  6  and  7* 

It  will,  be  seen  that  the  decay  times  of  the  two  modes  are  equal  when 
a  54O.  In  this  case  the  damping  ratio  of  the  oscillatory  mode  is  about 
O.58,  v/hich  is  satisfactory. 

The  value  of  w  was  chosen  to  be  equal  to  w0  in  this  analysis,  which 
is  equivalent  to  an  initial  demand  for  jet  angle  from  the  guidance  system  just 
equal  to  the  limiting  value,  The  effect  of  varying  the  natural  frequency, 

w,  of  the  autopilot  loop  is  considered  in  seotion  11. 

8  The  simulator  arrangement 

Initially  a  lab  oratory-built  simulator  was  used,  but  at  a  later  stage  the 
work  was  transferred  to  a  G-PAC  simulator.  The  equations  set  up  on  the  simu¬ 
lator  were  equations  (ll),  (12),  (13),  (l6)  and  (17),  the  shaping  network  in 
equation  (16)  being  a  simple  phase  advance :-* 


a  1  rad/sec,  w  a  1  rad/ 
results  are  shown  in 
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s  1  +  P®A 
Yl(p)  3  1  +  apTA  * 

A  block  schematic  of  the  simulator  arrangement  is  shown  in  Pig. 8. 

It  was  soon  found  that,  contrary  to  normal  experience,  the  presence  of 
phase  advance  had  a  destabilizing  effect  on  the  system  for  the  values  of  gain 
used.  The  phase  advance  network  was  therefore  by-passed,  and  with  the  para¬ 
meter  values 

wQ  =  1  rad/sec;  T^  s  2000  seo;  Tq  =  q:  sec; 

K,  =  1;  14,  =  1-6;  ^  =  550;'  I^A  »  £ 

-  */ 

the  simulator  records,  shown  in  Pig* 9  were  obtained  in  response  to  an  initial 
step  in  ^  of  magnitude  0.003  rad. 

These  results  confirm  the  rough  paper  analysis  given  above,  and  show  that 
the  system  can  be  satisfactorily  damped. 

Further  records,  all  in  response  to  a  0.003  rad  step  in  Ij,  showing  the 
effect  of  different  initial  conditions,  are  given  in  Pig. 10.  The  control  para¬ 
meters  have  the  same  values  as  given  above. 

In  Appendix  I  the  system  is  considered  from  the  frequency  response  point 
of  view,  and  an  analysis  of  the  effect  of  phase  advance  is  given. 

9  The  time  taken  to  reduce  an  initial  error 


In  the  ideal  system  (section  4  and  Fig. 5)  an  initial  flight  path  error 

is  reduced  linearly  with  time,  and  from  equation  (15)  it  is  seen  that  the 

time  required  to  reduce  to  zero  an  initial  error,  ,  is 

6 


Mi 


(24) 


In  practice  with,  T^  »  2000  seos  and  using  the  "optimum"  damping  arrange¬ 
ment  as  above,  it  was  found  from  the  simulator  work  that  the  time  required  is 
about  4  seos  more  than  that  given  by  equation  (24). 

The  simulator  results  are  compared  with  those  for  the  ideal  system  in 
Fig.11. 

10  Steady  state  errors 

The  use  ox  heading  angle  feedback  in  the  inner  loop  means  that  steady 
state  errors  oan  occur  due  to  (i)  changes  in  the  target  direction,  and 
(ix)  gyro  drift. 

Let  ijf1  be  the  measured  value  of  ljr  which  replaces  i|r  in  equations 
(16),  and  let  ijr1  -  »  A^r,  the  gyro  drift  angle.  Then  from  equations  (11 ), 

(l2),  (13)  and  (l6),  using  the  suffix  s  to  denote  steady  state  conditions, 

^  "  0 

*s  '  “  ^as 

*Ts  **  a  +SS  “  *is 

**  *  ~  *1  +  Vs  * 
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It  follows  that  the  steady  state  flight  path  error  is 


*ss  =  5 qfc  (*Ds  +  A*s>  • 


Thus  if  the  gyro  axis  is  along  the  assumed  sight-line  at  launch, 
represents  the  ohange  in  the  sight-line  during  the  controlled  flight  and  A^ra 
represents  the  steady  gyro  drift.  In  theory  a  correction  coxild  be  made  for  the 
first  part  of  the  error,  but  in  practice  the  change  in  the  sight-line  is  likely 
to  be,  at  most,  a  few  minutes  of  arc,  whereas  the  gyro  error  may  be  a  few 
degrees. 

Also,  »  Kj  fil’d  so,  approximately, 


-  £s\ 


From  Pig. 2,  the  line  error  due  to  this  is 


AL  st  R- .  ijr 


* 


R  .  Aijr 


where  R  is  the  range  from  cut-off  to  impact. 

Por  example,  if  Kj  «  1 ,  Kj  a  500,  R  a  2000  n.  miles,  then  a  gyro  drift 
of  7°  would  cause  a  line  error,  AL,  of  about  ■§•  n.  mile. 

K3 

The  factor  r*  will  be  called  the  DRIPT  FACTOR. 

K1 

11  The  inner  loop  gain 

In  the  work  described  so  far  the  inner  loop  gain,  Kj,  has  been  arbitrarily 
fixed  so  that  w,  the  natural  angular  frequency  of  the  autopilot  loop,  was 
equal  to  wQ.  The  effeot  of  varying  Kz  can  be  studied  by  considering  the 
stability  cubic  equation  (23)  again.  If  we  assume,  as  in  section  ~l,  that  a 
suitable  criterion  for  the  choice  of  is  equal  distribution  of  the  total 
damping  between  the  two  modes,  i.e.  a  «  Y>  this  leads  to  a  determinate  value 
of  Kj.  The  equations  to  be  solved  are 


3a  a  1 . 6w  -  •=- 

X1 

2  2 
0+2a  a  w 


0a  B  “  •  w  . 

H  T,  0 

1  '' 

These  have  been  solved  for  various  values  of  w,  with  w0  a  1  rad/seo 
and  T^  =  2000  seos,  the  results  being  shown  in  Pigs.  12,  13  and  14* 

It  will  be  seen  that  the  available  damping  is  virtually  independent  of 
but  that,  as  w  is  increased,  the  gain  and  the  angular  frequency,  *0, 

of  the  quadratic  mode  increase,  but  the  damping  ratio,  ,  decreases.  The 

drift  factor  improves  as  w  increases. 
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The  effect  of  increased  inner  loop  gain  on  the  overall  system  was 
studied  on  the  simulator,  and  a  simulator  reoord  is  shown  in  Pig.  15  for  the 
case  K-j  =  3,  Kg  =  3«4^,  Kj  =  2,800.  It  will  be  seen  that  there  is  no 
improvement  on  the  previous  records,  so  that  the  only  advantage  of  a  higher 
gain  appears  to  be  some  improvement  in  the  drift  factor.  A  disadvantage  of 
a  high  gain  which  has  not  yet  been  studied  in  detail  is  the  effect  of  natural 
vibrations  cf  the  missile  on  the  measurement  of  ijr  and  ^r. 

12  Intermittent  guidance 

In  the  case  of  the  radar  guidance  system,  the  guidance  information  is 
not  available  continuously  but  only  at  specified  intervals  of  time.  The 
information  may  also  be  late  to  the  extent  of  one  time  interval  unless 
special  prediction  computation  is  done  by  the  ground  computer.  The  effect  of 
this  sampling  and  time  delay  has  been  studied  theoretically  and  on  the  simu¬ 
lator,  and  details  are  given  in  Appendioes  II,  III  and  IV.  In  this  part  of 
the  study  the  thrust-to-mass  ratio  of  the  missile  during  the  vernier  stage 
was  assumed  to  be  -g-g  rather  than  ?g,  the  former  being  a  more  recent  estimate 

P 

of  a  realistic  value  for  —  The  main  conclusion  is  that  the  proposed 

control  system  is  satisfactory  if  the  sampling  interval  and  time  delay  are 
less  than  about  1  sec.  This  appears  to'  be  well  within  current  estimates  of 
the  computer  time  delay. 

13  Results  and  conclusions 

(i)  It  has  been  shown  by  simulator  work  and  paper  analysis  that  a  simple 
control  system  involving  heading  and  heading  rate  feedback,  and  with  limits 
on  the  azimuth  angular  error  input,  can  be  used  for  the  control  in  azimuth 
during  the  vernier  stage.  A  relatively  high  gain  system  is  recommended  in 
order  to  make  full  use  of  the  available  correcting  thrust. 

(ii)  If  K^ ,  Kj,  K-2  are  coefficients  of  flight  path  error  ijfe, 

- ^  - J-  *  - the  recommended  values 

Kj  »  1  rad/rad,  where 


assumed 
oonstant  . 


taken  to  reduce  to 


if  T, 

*  T— 1  +  (St) 

*max 


nesaxog  angle  y,  ana  neaaxng  rate  respectively, 
are  Kj  =0.27  £)  rads/rad.  Kg  =  i  •  °/wQ  rads/ rad,/ sec. 


missile  momentum 
thrust 


(in  sees) 


w_ 


thrust  x  moment  arm  / .  -2  \ 

moment  of  inertia  s  0  ) 


(iii)  With  the  above  control  coefficients  the  time 
zero  an  initial  error,  y_,  in  the  flight  path  is 


where  is  the  maximum  allowable  heading  angle  and  6t  is  approximately 

4  secs  when  T-j  =  2000  seos  (see  Pig.  11), 

(iv)  The  use  of  the  heading  term  in  the  control  equation  leads  to  steady 
state  flight  path  errors,  and  hence  to  a  resultant  line  error,  due  to  gyro 
drift.  'The  resultant  line  error,  AL,  due  to  a  gyro  drift  of  Afs  is 
given  by 


Hi  * 

-  12 


R  .  A  if 
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where  R  is  the  range  from  oat-off  to  impact. 

Thus  if  s  2000  secs  and  R  =  2000  n.  miles,  a  gyro  drift  of  7° 
causes  a  line  error  of  about  ^  n.  mile. 

(v)  Ream  (ii),  (iii)  and  (iv),  it  is  seen  that  if  the  thrust-to-mass 
ratio  is  reduced,  the  time  to  oorrect  an,  initial  error  will  be  longer,  but  the 
line  error  due  to  gyro  drift  will  be  less. 

(vi)  The  effect  of  sampling  the  flight  path  error  at  intervals  of  t  secs 
and  the  effect  of  a  finite  delay,  t,  in  the  correcting  signal,  have  been 
studied  for  the  case  T-j  a  l+OOO  seos,  w0  =  -g-  rad/sec  ( corresponding  to  the 
latest  estimate  of  thrust-to-mass  ratio  of  g/8,  as  mentioned  in  section  12). 
Provided  the  gain  is  slightly  reduced  below  the  value  given  by  (ii)  above, 
the  sampling  and  time  delay  have  virtually  no  eff eot  up  to  1  =  1  Beo. 

14  Further  work 

(i)  It  is  proposed  to  extend  the  investigation  of  azimuth  control  to  the 
main  stage,  where  the  variation  of  mass,  speed,  etc.  oannot  be 
neglected. 

(ii)  The  effect  of  (a)  flexural  oscillations  and  (b)  fuel  sloshing  on  the 
control  system  will  have  to  be  studied. 

(iii)  Work  is  to  be  done  on  -the  effects  on  the  control  system  of  various 

kinds  of  "noise",  and  on  Ihe  data  smoothing  required  to  minimise  such 
effects. 
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APPENDIX  I 


Prequenoy  response  analysis 


Taking  the  Laplaoe  transforms  of  equations  (ll),  (12),  (13)  and  (l4)> 


P  t 


ptm  =  ^  (5  +  2  -  im) 


2v 

~  wo  S 


(1+T^)^  =  ? 


2  =  -K^Cp)^*  (KgP  +  K3)i 


(28) 


Eliminating  f  from  the  first  two  equations 

7  2  2 
if  p  ~w 

jm  ^  *  o 

s  P2(1+T1p) 

Henoe,  substituting  for  ^  in  terms  of  ijrm  and  for  f'e  in  terms 

of  |E  in  the  final  equation  the  open-loop  transfer  function  is  obtained, 
namely 


i  a  -  yi(p)(p2-^02) 

*e  (l+T0P)(ltT1p)(p2+wo2K2p+wQ2K3) 

=  \(p),  say. 

With  the  figures  giver  in  the  text,  namely 

wq  =  1,  T1  .  2000,  Tq  =  i,  £,  =  1.6,  K3  «  1, 

and  with  no  phase  advance,  i.e,  X^(p)  =  1,  this  gives 

i l(p) - &n±dl - 

(l+ip)(l+2000p)(p +1#6p+1 ) 

The  gain-  and  phase-log  frequency  curves  for  Xq,( ju)  have  been'  computed 
and  are  shown  in  Pig»l6. 

Prom  these  curves,  it  is  seen  that  there  is  a  phase  margin  of  about  60°, 
which  is  ample  for  stability  purposes. 

In  order  to  determine  the  effect  of  phase  advance  on  the  system,  the 
characteristic  polynomial  was  studied. 


K,  »  55 0 
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The  character  is  tio  ^  polynomial  is  the  denominator  in  the  closed-loop 
transfer  function 


*t.(p) 


Lv 

[i+y'CpJJ 


Putting 


^(p)  -  1  +  IAP  , 


in  equation  (2$),  the  characteristic  polynomial  is  found  to  be 


where 


Ap^  +  Bp^  +  Cp2  +  Dp  +  E 

A  a  T  T.  • 

o  1 

B  =  To  +  T1  *  "/Wl  -  K1TA 
0  =  1  +wo2K2(To+T1)  +  w^T^ 

D  =.  Wo%  +  Wo%(S04fl^  *  K^w/ 

E  a  Wq2(K3+K1 )  . 


The  Routh-Hurwitz  stability  criterion  demands  that,  for  stability 

B  >  0 


BO-DA  >  0 

and  (BO-DA)  D-B^E  >  0  . 


With  the  values  for  wQ,  T^ ,  To>  Kg,  K^  given  above  the  first  condition 
beoomes  KjT^  ■<  2800.  Thus,  for  a  given  Kj,  a  sufficient  increase  in  the 
phase  advance  time  oonstant  causes  instability. 

Por  K^  a  1000,  the  stability  condition  is  found  to  be  T^  <  2.0.  This 

system  was  tested  on  the  simulator  and  was  found  to  go  unstable  at  about 

T^  a  1.3.  This  discrepancy  is  due  to  the  phase  advance  impurity  which  was  not 

included  in  the  above  analysis. 
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APPENDIX  II 

Intermittent  guidance  data 


II. 1  Introduction 

It  has  been  shown  that  a  suitable  azimuth  control  system  can  be  designed 
for  the  vernier  stage,  provided  that  the  guidance  information  is  available 
continuously,  i.e.  regarding  the  overall  system  as  a  continuous  servo.  In 
the  radio  system,  however,  the  guidance  signal  is  computed  on  the  ground  by  a 
digital  computer  from  data  provided  by  a  doppler-radar  system.  This  means 
that  the  guidance  signal  is  not  a  continuous  function  of  time  but  may  take 
only  those  values  appropriate  to  each  new  set  of  input  data  to  the  computer. 
Suppose  that  there  is  an  interval  t  between  each  new  input  to  the  computer; 
then  guidance  signals  are  available  at  a  frequency  of  1/r  and  the  error 
quantity  may  only  be  measured  once  in  t  secs.  This  fact  means  that  the 
system  becomes,  in  effect,  a  sampling  servo,  and  there  is  no  guarantee  that 
the  suggested  oontrol  system  will  be  satisfactory  in  these  circumstances. 

A  second  effect  which  arises  from  the  use  of  a  digital  computer  is  the 
delay  introduced  between  measurement  of  the  data  and  transmission  of  the 
appropriate  signal  to  the  missile,'  in  effect  the  information  is  late  by  the 
amount  of  time  taken  to  perform  the  computation.  The  effect  of  this  delay 
has  been  examined  both  on  the  simulator  and  in  theory. 

II.  2  Intermittent  data  -  no  delay 

In  the  first  instance,  the  effect  on  the  missile  performance  of  sampling 
the  information  was  studied,  on  the  assumption  that  the  delay  due  to  computa¬ 
tion  is  negligible.  It  is  probable  that  the  data  rate  available  from  the 
computer  will  not  be  lower  than  1  per  second,  and  not  higher  than  5  per 
second,  so  that  the  simulator  investigation  consisted  of  examining  the  missile 
behav.iour  as  the  sampling  period,  t,  varied  between  y5  and  1.  It  is  neces¬ 
sary  to  maintain  the  loop  closed  between  sampling  instants  by  clamping  the 
guidance  signal  to  the  value  it  had  at  the  previous  sampling  instant.  Thus 
if  lJicCt)  is  the  output  of  the  clamp  and  ^r£(t)  is  the  guidanoe  signal, 

*Q(t)  *  ^(a*)  (a  =  0,1,2....) 

for  nr  <  t  <  (n+1  )-r  . 

The  simulation  with  intermittent  data  has  been  done  for  a  thrust-to-mass 
ratio  of  fg,  rather  than  the  used  hitherto  in  this  note,  in  line  with  more 
recent  thoughts  on  the  vernier  stage.  As  is  to  be  expected,  after  the  inner 
loop  has  been  appropriately  adjusted,  the  lower  available  control  acceleration 
roughly  doubles  the  time  taken  to  correct  an  initial  error.  Thus,  when  the 
error  information  is  continuous,  the  time  taken  to  correct  a  flight-path  error 
of  3  mils  is  approximately  30  secs  rather  than  16  secs  as  obtained  previously. 

The  same  simulator  arrangement  was  used  as  in  the  continuous  case,  but 
the  loop  was  broken  at  X  (Pig.8)  and  the  sampling  device  inserted.  The 
method  used  for  simulation  of  the  "sampling  and  clamping"  is  as  described 
in  Ref.  1 ;  further  details  are  given  in  Appendix  IV.  With  the  para¬ 
meters  chosen  to  give  satisfactory  performance  in  the  continuous  system, 
the  missile  response  to  a  3  nil  step  in  was  examined  for  *  =  1/5, 

t  =  2,  end  x  =  1.  It  was  found  that  the  intermittent  nature  of 
the  data  tends  to  destabilize  the  system  slightly,  but  even  when  t 
is  as  large  as  1  sec,  the  destabilizing  effect  is  small.  This  has 
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been  confirmed  by  paper  analysis  using  the  pulse  transfer  function  which 
shows  that  instability  is  delayed  until  ir  is  greater  than  about  1 6  seconds 
(Appendix  III). 

II.  3  Intermittent  data  with  delay 

The  above  simulation  does  not  take  acoount  of  the  delay  arising  in  the 
digital  computer  due  to  the  finite  time  required  for  computation,  and  the  simu¬ 
lator  was  therefore  modified  to  include  a  delay  of  length  v,  where  ^/x  is 
the  sampling  frequency.  The  method  employed  is  described  in  detail  in 
Appendix  IV. 

As  is  to  be  expected,  the  computing  delay  reduces  the  stability  of  the 
system.  For  <ir  »  1  sec  the  system  is  markedly  03oi.lla.toiy  but  still  stable; 
in  Appendix  III  it  is  shown  theoretically  that  instability  ooours  when  x  is 
greater  than  about  7  secs. 

In  an  attempt  to  regain  the  well-damped  behaviour  shown  in  the  case  of 
continuous  data,  the  gain  constant,  Kj,  was  varied.  It  was  found  that  for  each 
sampling  frequency,  i/v,  it  is  possible  to  choose  ah  "optimum"  value  for  Kj, 

-  the  word  "  optimum"  being  used  to  mean  that  for  an  initial  step  in  ^  the 
time  of  settling  is  a  minimum  -  and  that  the  value  of  the  "optimum"  Kj  for  a 
given  x  is  very  nearly  proportional  to,  x  (Fig.  17).  Thus  it  appears  that, 
provided  a  suitable  loop  gain  is  chosen,  the  performance  of  the  suggested  con¬ 
trol  system  is  only  slightly  degraded  when  the  error  information  is  sampled  and 
delayed.  Simulator  records  of  responses  to  a  3  mil  step  in  \jsji  are  given  in 
Figs.  18,  19. 

II. 4  Conclusions 

1.  Although  sampling  of  the  error  data  reduces  the  stability  of  the  system, 
this  reduction  is  very  small  for-  sampling  periods  of  up  to  1  sec. 

2.  The  effect  of  a  one  period  delay,  xf  in  the  main  feedback  loop  is  to 
further  reduce  the  stability;  this  effect  is  small  for  x  <  ■§•  sec,  but 
becomes  marked  at'  x  =  1  sec. 

3«  The  instability  caused  by  the  introduction  of  the  delay  can  be  compen¬ 
sated  for,  to  a  large  extent,  by  reducing  the  loop  gain,  and  the  necessary 
reduction  in  loop  gain  is  approximately  proportional  to  the  sampling 
period,  i.e,  to  the  delay. 

4.  The  system  suggested  for  azimuth  control  during  the  vernier  stage  is 
satisfactory  even  if  the  error  data  is  sampled  at  a  frequency  of  ^/x 
samples/sec  and  delayed  by  x  secs,  provided  that  x  is  not  more  -than 
about  1  second. 
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APPENDIX  III 


Pulse  transform  analysis  of  system  with  intermittent  data 


(i)  No  delay 

Following  the  method  of  Ref.1  we  can  split  up  the  forward  path  of  the 
loop  into  two  filters,  such  that  information  passes  from  the  first  to  the 
second  filter  only  in  samples. 

We  split  up  the  clamp  plus  missile  system  into  two  components,  the 
first  filter  representing  an  imaginary  sampling  device,  and  the  second  having 
the  Laplace  transform  p~'l  Yp(p),  where  Ip(p)  is  given  by  equation  (29). 
[The  justification  for  this  step  is  given  in  the  first  appendix  of  ref.I.J 


Then 


*l(p> 


V^q2-p2) 

(1+^p)  (p2+wo2K2p+wo2K3  ) 


and  so 


p"1  tt.(p)  = 


xr  /  2  2\ 

Vwo  “P  > 


p(l+T1P)(p2+wo2X2p+wo2K'3) 


-P.».09J.Z5...,1.0.r.187^. 


+  1.1 20p3  +  0.5003P2  -I-  0.000125P 


This  has  the  weighting  function: - 

W(t)  =  [p"1  IL(p)} 


i.  e. 


W(t) 


750 


D 


-  0.99 


t 

4000 


+  e 


-0.56t 


(0.416  cos  0.42t  -  O.316  sin 


0.42t) 


Thus,  the  weighting  sequence  is:- 


W(nv) 


750 


1  -  O.99  e 


nv 

4000 


-0«56nv 

+  © 


(0.416  00s  0.42nr  -  0.316  sin  0»42mr) 


Nov/  let  the  pulse  transfer  functions  of  the  two  filters  be  u(z)  and 
V(z)  respectively,  so  that  the  total  loop  transfer  function,  W(z),  is 
given  by 
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W(z)  =  U(z)  .  V(z)  . 


Then 


r — i  f — i  /.  2*000 \n 

V(z)  »  750  \  z  n  -  750  \  (2—z — ■) 


n=o 


nao 


+  V"'  fe“°^6nt(0<41g  QQg  Q.42nT  -  0.316  sin  0»2*2nT  J 


nao 


Performing  the  summations,  we  obtain:- 
V(z)  a 

v  z-1  za-1 


where 


O.M6(z3r  -  zB( 0.4-1  go  *  0«3l6s> 

(z(3)2  -  2zj3C  +  1 


a  a  e2*00 
0  a  e0'56* 

0  a  OOS  0.2*2* 

S  a  sin  0.2*2*  . 


Then 


w(z)  a  u(z)  .  v(z) 

=  (W*1)  •  V(z) 


So 

w(z)  «  252£a=dl  +  (zm1)  -  P(o.fci6o  ±  0.31 6s  n  # 

L  (zp)2  -  2zpC  + 1  J 

Now  for  stability  of  the  closed-loop  system,  all  the  roots  of  the 
characteristic  equation  of  the  overall  transfer  function  must  lie  inside  the 
unit  circle  (Ref.i  -  Appendix  i). 


is:- 


The  characteristic  equation  for  the  overall  transfer  function, 
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z^(a(32  +  0.4aj32) 

+  z2(~(32  -  2a(3C  +  750(a-l)p2  -  0.4p2(l+a)  -  a(3(0.4l6c  +  0.316s)) 
+  z(<x  +  2pC  -  750(0-1)  2|3C  +  O.k-f  +  (1+a)  p(0.4l6c  +  O.316S)) 

+  (75°(a-1 )  -  1  -  (3(0. 416C  +  O.316S))  =.  0  . 


It  is  easily  shown  that  to  determine  whether  the  roots  of  this  equation 
lie  inside  the  unit  circle  it  is  only  neoessary  to  transform  the  interior  of 
the  unit  circle  in  the  z-plane  into  the  left  half  of  the  p-plane,  using  the 
transformation  (for  example):-* 


Then  the  Routh-Hurwitz  criterion  for  roots  in  the  left  half  plane  oan  be 
applied  to  the  transformed  equation. 


As  an  approximation  we  set  a  =  1,  exoept  in  terms  containing  750  (a~1 ) ; 
then  the  coefficients  of  the  transformed  equation  beoome:- 


Aq  =  750 (a-1 )  [P2  -  20(3  +  1] 

A1  =  [2  +  750(a-l)]p2  -  [40  -  1500(cc~l)C]p  +  [2  -  2250(a~1 )] 

■^2  “  [5.6  -  750(a-l)]p2  -  [4Y  -  1500(a-1  )C](3  -  [4  -  2250( a-1 ) ] 

A3  «  [3.6  -  750(a-*1  )]02  +  [4(r+C)  -  1500(a-1  )]p  +  [2  -  750 (a-1 )] 


where  Y 


0.4160  +  O.316S. 


Substitution  of  trial  values  for  v  leads  to  the  condition  for  stability 
that:- 


v  <  ~  16  secs  for  stability. 


(ii)  Delay  in  forward  path 

It  oan  be  shown  that  when  there  is  a  delay  of  m  periods  in  the  error 
path,  the  overall  pulse  transform  becomes 


z  m  .  W(z) 

1  +  a"®  .  W(z) 


Thus  for  a  delay  of  one  period,  the  characteristic  equation  becomes 


z  +  W(z)  e  0 


and  again  the  roots  mist  lie  within  the  unit  circle  for  stability. 
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In  detail,  the  equation  becomes:- 
z\af]  +  z3  [-p2  -  2a(3C  +  O.4a023 

+  z2  [a  +  2|30  +  750(l->a)p2  -  0,4p2(l+a)  -  af3(0.4l6C  +  0.316s)] 
+  z  [~  1  -  750(a-l)  20C  +  0.4P2  +  (1+a)j3  (0.4160  +  O.3I6S)] 

+  [750(a-l)  -  p( 0.4160  +  0.316S)]  «  0  . 

Evaluation  of  the  coefficients  of  the  ooixesponding  transformed  equation, 
and  substitution  of  trial  values  of  t  leads  to  the  condition: - 

x  <  ~  7  seos  for  stability. 
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APPMDJX  IV 
Simulator  details 


IV.  1  General  $ 

Most  of  the  simulation  was  performed  on  a  single-cabinet  G-PAC  computer, 
consisting  of  20  drift-corrected  D.C.  amplifiers  and  built-in  power  supplies. 
The  computing  networks  are  mounted  on  plates  which  plug  in  to  the  amplifiers, 
so  that  computing  impedances  are  readily  changed  if  neoessazy.  Interconnec¬ 
tions  are  made  from  the  fronts  of  the  oomputing  plates  and  signal  bus -bars 
are  available  at  the  side  of  the  cabinet,  as  are  the  computing  voltages 
(±  100V)  and  step  voltages  (±  100  V).  Pacilities  are  available  for  intro¬ 
ducing  initial  conditions  onto  the  amplifiers  when  used  as  integrators. 

IV. 2  The  sampling  device 

The  simulation  of  the  sampling  device  is  basioally  that  used  in  Ref.  1 . 

An  asymmetric  multivibrator  drives  a  high  speed  relay  whose  moving  contact 
carries  the  error  signal  (Pig.  20 );  for  most  of  the  multivibrator  period  (the 
"space")  the  moving  arm  makes  with  an  isolated  fixed  contaot,  but  during  the 
relatively  short  "mark"  period  contact  is  made  with  the  other  fixed  contaot 
which  is  connected  via  a  small  resistor  to  the  storage  capacitor  (clamp). 

(The  purpose  of  the  small  resistor  is  to  ensure  that  the  cathode  follower  of 
the  previous  amplifier  is  not  called  upon  to  supply  too  much  current). 

Since  the  relay  makes  contact  for  a  finite  time  -  approximately  the 
"mark"  period  of  the  multivibrator  -  the  arrangement  is  not  a  true  sampling 
device,  i.e.  its  output  is  not  in  the  form  of  6-f 'motions.  However,  provided 
that  the  mark-space  ratio  of  the  multivibrator  is  less  than,  say,  0.Q5,  a 
good  approximation  to  a  sampling  device  is  obtained. 

IV. 3  The  delay  mechanism 

The  production  of  a  delay  equal  to  one  sampling  period  is  obtained  by 
using  an  extension  of  the  sampling  simulator  previously  described.  An 
asymmetric  multivibrator  drives  a  Type  3000  heavy-duty  relay;  during  the 
short  "mark"  period  this  belay  closes  a  circuit  putting  24V  on  the  coil  of  a 
standard  four-bank  double-arm  uniselector.  The  input  voltage  (error  signal) 
is  oonneoted  to  the  first  wiper  and  the  output  lead  is  connected  to  the 
second  wiper;  a  storage  capacitor  is  connected  (through  a  small  resistor) 
to  the  first  contact  of  bank  1  and  the  output  of  this  condenser  is  taken  to 
the  third  contact  of  bank  2;  similarly  a  second  storage  condenser  is  con¬ 
nected  between  contact  2  of  bank  1  and  contact  4  of  bank  2,  and  so  on.  In 
this  case  five  storage  condensers  were  used,  and  every  fifth  contact  on  each 
bank  was  comaoned.  The-  mode  of  operation  is  as  follows.  Suppose  that  all 
the  capacitors  have  no  charge  on  them  and  that  at  t  =  0  the  wipers  are  on 
contact  1;  then  at  t  n  r  the  wipers  move  on  to  contact  2  and  the  error 
signal  is  applied  to  contact  2,  bank  1.  The  output  (i.e.  the  voltage  on 
contact  2,  bank  2)  is  zero,  since  this  is  the  voltage  on  the  capacitor  whose 
input  is  contact  50,  bank  1  (3  contact  5*  bank  l).  Now  at  t  =  2-v  the 
wipers  move  on  to  contaot  3  and  the  voltage  left  on  contaot  2,  bank  1,  is 
that  appropriate  to  the  error  signal,  at  t  »  and  not  t  »  t.  The 
voltage  picked  off  oontact  3>  bank  2,  is  the  output  of  the  capacitor  whose 
input  is  contact  1 ,  bank  1 ;  i.  e.  the  output  at  t  a  2t  is  the  error  voltage 
appropriate  to  tat.  Then  at  t  n  Jt  the  wipers  move  to  contact  4  and 
in  this  case  the  output  wiper  picks  off  the  voltage  frem  the  capacitor  whose 
input  is  contaot  2,  bank  1.  This  voltage  i3  that  corresponding  to  the  error 
at  t  =  2r  (see  above)  and  so  a  delay  between  input  and  output,  of  length 
T,  is  obtained. 
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In  final  form,  the  sampling  unit  arid  delaying. unit  are  ootibined  into  one 
unit,  and  a  switching  arrangement  is  provided  so  that  the  output  of  the  unit 
may  be  sampled  but  not  delayed,  or  sampled  and  delayed  by  one  sample  period. 
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FIG.  I  &  2. 


MISSILE 


FIG. 2  RELATION  BETWEEN  LINE  AND  ANGULAR 

ERRORS. 


FIG.3.  BLOCK  DIAGRAM  OF  COMPLETE  SYSTEM. 
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FIG.6&7. 
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FIG.  8. 


FIG. 8.  BLOCK  SCHEMATIC  OF  SIMULATOR  ARRANGEMENT. 
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FIG.9(o-<0 


(b)  HEADING  ANGLE. 


(C)  JET  DEFLECTION. 


(d)  FLIGHT- PATH  ANGULAR  RATE. 

FIG.9.  (a-d)  SIMULATOR  RECORDS  -  RESPONSES 
FOR  STEP  DEMAND  OF  TJfT  (003  RAD.)  K ,  -  550, 
K2 *  1-6,  Kj-I,  NO  PHASE  ADVANCE  f0  =  %*O. 
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FIG.II.  TIME  TAKEN  TO  REDUCE  INITIAL  ERROR, 

PLOTTED. 

(q)  as  a  function  of  Tjrt  (with  i|rMAy  =  o-s  rad.) 

(b)  AS  A  FUNCTION  OF  ljfMAX  (WITH  =  3  MILS.) 

(c)  CURVE  (a)  FOR  THE  I0EAL  SYSTEM  'l 

(d)  CURVE  (6)  FOR  THE  IDEAL  SYSTEM  J  WITH  T'  =  2000  SECS. 
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FIG  12.  EFFECT  OF  INNER  LOOP  GAIN  ON  GAIN 
REQUIRED  TO  GIVE  EVEN  DAMPING  DISTRIBUTION. 


FIG  13.  EFFECT  OF  INNER  LOOP  GAIN  ON  SYSTEM 
FREQUENCY  AND  DAMPING  RATIO  WITH  THE 
GAIN  CHOSEN  TO  GIVE  EVEN  DISTRIBUTION 

OF  DAMPING. 
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FIG.  14.  EFFECT  OF  INNER  LOOP  GAIN  ON 
THE  DRIFT  FACTOR  K,/K3  WHEN  THE  DAMPING 
IS  EVENLY  DISTRIBUTED  BETWEEN  THE 

MODES. 
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FIG.l5(q*b)  RESPONSES  WITH  A  HIGHER 
INNER  LOOP  GAIN. 

K,  =  2,800,  Ks=  3*46,  K3=3. 

NO  PHASE  ADVANCE. 
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FIG.I7 


FIG.  17  VARIATION  OF  OPTIMUM  GAIN  WITH 
SAMPLING  INTERVAL  AND  DELAY  TIME  OF  T. 
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(g)  X  *  I  SEC.-  DELAYED  DATA -OPTIMUM"  K, 
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FIG.  19  (e-h)  COMPARISON  OF  HEADING  RECORDS 
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FIG. 20.  CIRCUIT  OF  SAMPLING  DEVICE. 
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